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p Abstract 

■ The restoration of chiral symmetry at finite density and/or temperature is 

investigated in a wide class of one-gluon exchange models in the instantaneous 
approximation. If the effective quark interaction is less divergent than 1/k 2 
for small momentum transfer k, we obtain Gaussian critical exponents for 
the chiral phase transitions at finite temperature and density, respectively. In 
the opposite case, for an interaction diverging faster than 1/k 2 in the infrared 
region, a qualitative different behavior of the quark self-energy near the critical 
Fermi momentum k c and the critical temperature T c , respectively, is observed. 
In the first scenario, we find k c ~ 2 In 2 T c , which compares well with recent 
data from QCD lattice simulations. 
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1. Introduction 

The properties of hadrons in a hot and/or dense medium are of general interest, 
since they can be experimentally tested in heavy ion collisions. These experiments 
might provide informations of the widely unknown low energy sector of QCD. Fur- 
thermore, a detailed understanding of the quark ground state at finite density is 
necessary for a wide range of applications, e.g. to describe compact star matter ||T|. 
In these applications, the restoration of the spontaneously broken chiral symmetry 
is of particular interest. For medium values of the density compared with the criti- 
cal density, a powerful argument was developed which predicts a scaling of certain 
hadronic observables with the density 0. This so-called BR-scaling yields phe- 
nomenological reasonable results for densities up to the critical density. Close to the 
critical values, physics is described by another approach; it was argued long ago 0, 
that the chiral transition at finite temperature is governed by dimensional reduc- 
tion and universality. In particular, four-dimensional QCD should belong to the 
universality class of the three-dimensional a-model with the corresponding global 
symmetries 0. Recently, Kocic and Kogut argued, that the presuppositions which 
support the universality argument, are not satisfied. In particular, they observe that 
the critical exponents of the chiral transition in the three-dimensional Gross-Neveu 
model are of Gaussian type, instead of those given by the Ising model, as predicted 
by the universality argument [f|]. The chiral transition at finite density is less known, 
which is partly due to conceptual and practical difficulties of the numerical simula- 
tion at finite densities |5[]. It is tacitly assumed that the critical value and indices 
of the chiral transition at finite density resemble those of the transition at finite 
temperature. 

In this letter, we study the chiral transition at finite density and/or temperature 
in a wide class of constituent quark models with effective one-gluon exchange. We 
will find that Gaussian critical exponents are the generic case, and will discuss 
the circumstances which imply a correspondence of the chiral transition at finite 
temperature and finite density. 

2. Model description 

Starting point of the one-gluon exchange models is the QCD Dyson-Schwinger equa- 
tion for the quark propagator S(p), which in Landau gauge and in Euclidean space- 
time reads 



if 1 + m + 



/ 



d A k 



r«( P ,p-k)s(k) lu t a D; b v (p-k) = s-\ P ), 



(i) 



(27T) 4 



where the Euclidean quark propagator is of the form 
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Equation ([!]) must be supplemented with Ansatze for the full gluon propagator 
D"^(p—k) and the quark-gluon vertex function T^(p, p—k). The variety of models 
|7|, ||, ||, ITOf-and this is not a complete list-differ in the particular choice of D°^ v and 
Y a . These choices must match the well known behavior of the perturbative Green's 
functions at large momentum transfer, i.e. 



D tM = ..L is , for # » Kl CD . (3) 



167T 2 1 / _k t k, 

lllnF/A 2 ,^ k 2 { P 
Sometimes, the low-energy behavior of the gluon-propagator is chosen as 



KM cx 8 ab i- A ( 6,., - ^ ] for k'^0. (4) 




This choice produces a linear confining potential between static quarks. It was 
observed, that the ansatz (f|) yields a model which is one-loop renormalizable and 
realizes confinement by infrared- slavery [H], which enforces the quark propagator 
to vanish. In any case, we will assume that the interaction under consideration is 
strong enough to break chiral symmetry spontaneously. 

A different confining mechanism arises when one assumes that the Yang-Mills vac- 



uum is realized by random color electro-magnetic fields pT| , which induce a complex 
constituent quark mass. In this case, the quark propagator is finite, but quarks do 
not exist as asymptotic states, and furthermore the Landau pole of the correspond- 
ing perturbative gluon-propagator is screened due to the presence of a dynamical 
mass |l2j . This implies that the behavior 

« ^ (V - ^) for fc2 -°- ( 5 ) 
is consistent with phenomenology. 

In this letter, we need not to specify the full gluon-propagator beyond the constraint 
(|]), but will refer to the different specific scenarios Q) and ([5]) when needed. 
In order to stress the generality of our main result, we investigate to different Ansatze 
for the vertex function r?., i.e. the bare vertex 

r>,g) = (6) 

which is the simplest choice, and a more sophisticated vertex, i.e. 

r^p,*) = t^T.ip^k^ + ip. + k^ + miP 2 ^ 2 ) (7) 
+ (^ + ^)(^-^)r 3 (^A; 2 )-^ + A; M )r 4 (^P)[s(fc 2 )-S(j9 2 )" , 



which is motivated from Taylor-Slavnov identities P, 13]. The last term on the 



right hand side of (0), which is proportional to the self-energy, ensures that the 
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Ansatz does not explicitly break chiral symmetry. From dimensional arguments, 
one observes that the functions T 2 ..a decrease like 1/p 2 and 1/k 2 for large values for 
p 2 and k 2 , respectively. 

Using the Ansatz (^), the Dyson- Schwinger equation decays into two equations for 
Z(p 2 ) and the quark self-energy S(p 2 ). The equation for the latter can be cast into 
the form 

d 4 k E(/c 2 ) 



5V 



m + 



■V(k,p) 



(2tt) 4 Z 2 (k 2 )k 2 + Y?{k 2 ) 

the functions H and T> are 
1 . 



H(p 2 ) 

where for the case of the bare vertex 

H{ P 2 ) : 

V{p,k) = ±D%(p-k), 

while in the case of the more complicated vertex (|7|), these functions are 

d 4 k Z{k 2 ^ 



H(p 2 ) 
V{p, k) 



1 + 



(2vr) 4 Z 2 (k 2 )k 2 + T 2 (k 2 



(9) 
(10) 

(11) 
(12) 

(13) 

Note that the function H(p 2 ) in (|ll|) is ultra-violet finite, since T>i(p,k) decreases 
like 1/k 4 for large k. 

In order to address the renormalization of the Dyson-equation, we exploit the fact 
that the self-energy S(/c 2 ) and the gluon induced interaction T> match the pertur- 
bative behavior at large momentum transfer, i.e. 

^. l9 . const. m const. , , 

nk 2 ) « 7 I ^ . V « TTT-TTTT^- , (14) 



where 



Z(fc 2 )T 4 (p 2 , k 2 ) + 2T 2 (p 2 , A: 2 )j P/iJ D> + -T^p 2 , k 2 )D 



V^k) = Up 2 ,k 2 )pMp-k)p 



In k 2 / h.Q CD 



k 2 In k 2 jk 2 QCD ' 



where c? m > is the anomalous dimension of the quark mass. The generic behavior of 
Z(k 2 ) turns out to be weakly momentum dependent and asymptotically approaches 
a constant for large k 2 . In fact, Z(u 2 ) is renormalized to 1 at some large subtrac- 
tion point [x in perturbative QCD with mass shell renormalization |TJ]]. This large 
momentum behavior of Z(k 2 ) and S(/c 2 ) implies that the momentum integration in 
the Dyson-equation (§|) is ultra-violate finite. If one nevertheless cuts off this mo- 
mentum integration by an 0(4)-invariant cutoff A, the behavior of a self-consistent 
solution of (H) enforces the bare quark mass m to vanish for A — > oo as (see e.g. 

0) 



m(A 2 



1 



[lnA 2 //i 2 



m R (a) 



(15) 
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where /i is the renormalization point and wir is the renormalized quark current mass. 
Here we are only interested in the chiral limit implying 

^R^f) = — > m = (16) 

from the very beginning. 

In the following, we will adopt further simplifications which we believe not to change 
our conclusions. We will assume that the function Z(k 2 ) is weakly dependent on 
the momentum k 2 for the models which we consider here, implying that we may 
safely approximate Z(k 2 ) = 1. In fact, one realizes that Z(k 2 ) = 1 is a solution 
of the Dyson-equation, if one neglects vector-contributions to the Dyson-Schwinger 
equation after a Fierz rearrangement of the interaction. 

Finally, we assume a separation of the fermionic and gluonic energy scales Ef and E g , 



i.e. Eg S> Ef. This assumption is motivated by lattice results ||16|| , which show that 
the lowest-lying excitations of the Yang-Mills vacuum, the glue balls, have energies 
larger than w 1.5 GeV, while the energy scale of the lowest quark modes is set by the 
constituent quark mass of ~ 300 MeV. The above presupposition has two immediate 
consequences; firstly, the contribution of the temperature dependent gluon exchange 
to the loop integration in the Dyson-Schwinger equation is of order exp{— E g /T}, 
which is small compared with the fermionic part exp{— Ef/T}. This implies that 
we may neglect the temperature dependence of the interaction. In fact, this idea 
is supported by simulations of lattice QCD, where one observes that the gluon 
condensate is weakly temperature dependent throughout the chiral transition || 
suggesting that the chiral transition is solely driven by the temperature dependence 
of the quark loop. Secondly, the dominant contribution of the ko integration in the 
loop integral of the Dyson-Schwinger equation stems from the region k Q « Ef^\. This 
implies that we may neglect the ko dependence in the gluon exchange, and we end 
up with an instantaneous approximation. 

3. Restoration of chiral symmetry at finite density 

In order to investigate density effects, we introduce a chemical potential /i, which 
enforces a non- vanishing vacuum density p = (Q\q^q\Q). In instantaneous approxi- 
mation, the quark self-energy £ only depends on the three momentum k, implying 
that the fco-fntegration in the Dyson-Schwinger equation (|§D and in H(p 2 ) (|ll]) can 
be done, i.e. 

Hfip 2 ) h*>k) (2tt) 3 2^]k 2 + Y 2 (k 2 ) 

r d 3 k 1 -> 

H f m = i+ 1 * 3 ,= _ pi(p,ao, (is) 

Jk2 > k2 f ( 2?r ) 2^Jk 2 + T 2 {k 2 ) 



1 This argument becomes obvious, if one performs this integration in the complex plane by 
Cauchy's theorem. 
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where the subscript "p" indicates the finite density case. Furthermore, the Fermi 
momentum kf is related to the density by k 3 = ir 2 p for N c = 3 colors and Nf = 1 
flavor. In the case of the bare vertex @, the function Hf is again Hfip 2 ) = 1 and 
thus independent of the Fermi momentum. 

We first derive an equation for the critical Fermi momentum k c , at which chiral 
symmetry is restored. To this aim, we apply the so-called bifurcation method fIT|| 
and exploit the fact that at the critical density, a second solution to the Dyson- 
Schwinger equation ( |17|) exists besides the trivial one E = 0. At the critical density, 
the second solution differs from the trivial one only by an tiny amount implying 
that its momentum dependence /(jo 2 ), where Efj) 2 ) oc with /(0) = 1, can be 

calculated from the linearized equation, i.e. 

f(?) = TFT^rT / „ 7^ A=V&P) f(H (19) 



Hoip 2 ) h*>k* (2tt)3 2 ^ 
where 

r ri 3 k 1 

H " (f) - 1 + L } m^ VMk) ■ (20) 

This equation is essentially the Bethe-Salpeter equation describing fluctuations 
around the trivial vacuum at the critical point. Since the linearized equation is 
homogeneous in E, it does not fix the magnitude of the self-energy. The critical 
density is obtained from (|T9| ) by demanding that f{p 2 ) is an eigenvector of the 
integral- kernel on the right hand side of (^) with eigenvalue 1. 
In order to fix the size of the self-energy for a Fermi-momentum close to k c , we study 
the next to leading order of the small amplitude expansion of the Dyson-Schwinger 
equation (|Hj) 



i.e. 



1 f d 3 k E(£ 2 ) 

^ = L',wfi/i mk) (2l) 

1 f d 3 k S(fc 2 ) g f d 3 q S 2 (^) -> 
+ 2Hl(f) J*Z* f W 2 Jp [P ' } U>*, W 2(^)3/2 V W> k > 

1 /■ d 3 fc E 3 (£ 2 ) . w A 



2H (p 2 ) Jkz>kj (2tt) 3 2(P)3/2 

We will solve this integral equation approximately. On first realizes that for large 
external momentum p 2 the last two terms of (|2l|) are irrelevant, since the functions 
Dn ) (p, k) are strongly peaked for k 2 m p 2 implying that the these last two terms are 
suppressed by a factor E 2 jp 2 compared with the first term in (JUT). We therefore 
conclude that for p 2 ^> k 2 , the momentum dependence of the solution Efjo 2 ) of 
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is given by f{p 2 ) from (0). This motivates the approximation 



(22) 



where Em must be calculated from (21) . The procedure implies that the ansatz 
(B^) is a good approximation to the full solution of (EHI) at least at small and at 
large momentum p 2 . Using the implicit expression ([19|) for k c , one obtains to leading 



order in k 2 c — k 2 from 



J (i) 



d 3 k 



11) at p l = 

rf 3 g V{0,k)f 3 (k 2 ) 



fc2>fc| (2vr) 3 y g 2> fe 2 (2tt) 3 2(P) 3 / 2 



(23) 



S(k-q)-e 



fc2^ 3/2 p,i()./.'j 

9 s 



fc 2 <* 2 _fk_ f(k 2 ) 
k 2 >k* (27r) 3 



P(0, k) 



where e = for case of the bare vertex ([]) and e = 1 for the case of the improved 
vertex (|7|). If one assumes that the function T^p 2 ^ 2 ) is not singular for p 2 — > 
(which is the generic behavior |13|), the function T>i(p 2 = 0,k 2 ) vanishes implying 
that our result ( p3"D is widely independent of the ansatz for the vertex function. The 
only dependence on the particular choice enters via the function Ti(p 2 , k 2 ) through 
the projected gluon propagator V(p, k). 
The integral at the right hand side of (|23|) is proportional to 



duf(u) -T x (0, u)D™(u). 



(24) 



We therefore obtain finally 



E (i) 



klf{kl)T x {Q,kl)D™{k 



An 2 J, 



d 3 k T!(0,fc 2 )D-(fc 2 )/ 3 (fc2) 



fc 2 >fc 2 



(2tt)3 



( fc 2)3/2 



1/2 



\ 



k 2 

1 ^> 

k 2 



(25) 



which corresponds to a Gaussian critical exponent 1/2. 



4- Restoration at finite temperature 

The quark interaction of the one-gluon exchange model is motivated from the glu- 
onic interaction of QCD and therefore temperature dependent. However, numerical 
simulations of lattice QCD indicate that the gluonic properties only vary weakly, 



if the quark sector passes the chiral transition |18[. We therefore neglect here any 
explicit temperature dependence of the interaction in (P, but study the tempera- 
tures effects in the quark sector, since those apparently play the dominant role at 
the chiral transition. 



7 



The standard method to calculate the impact of temperature on Green's function 
is the imaginary-time formalism |B|, i.e. one confines the configuration space of 



the fermionic fields to the configurations which are anti-periodic in the Euclidean 
time direction with a periodic length 1/T with T being the temperature. At finite 
temperature, the integration over the zeroth component of the Euclidean momentum 
in the terms of (|8]) is replaced by a discrete sum over Matsubara frequencies. In the 
case of the instantaneous approximation, this sum can be explicitly evaluated and 
yields 



E(p 2 ) = — — - — - tanh — v J V(k,p) , (26) 

KF J H T (f)J (2vr) 3 \2TJ 2E{k 2 ) 1 ' 1 ' 

H T (f) = 1 + e / tanh ( -=-) ^V±{p,k) , (27) 

^ ; J (2vr) 3 V2T7 2E(k 2 ) 



where the subscript "T" indicates the finite temperature case, where E(k 



2^ 



k 2 + Ti(k 2 ) and again e = for the case of the bare vertex (|6|) and e = 1 oth- 
erwise. The close correspondence of the density and temperature dependence of the 
model stems from the similarity of the equations ( p6| - p7| ) and flT7j-|T8|) respectively; 
whereas the tanh in (^6||27]) rapidly approaches 1 for E ^> T, small momenta at 
E CT are suppressed. The same behavior is observed in the kernel of the integral 
equations (|17| - |i~8|) in the case of finite density. The only difference is that in this 
case the screening of low momenta is done by a step function due to Pauli blocking, 
whereas low momenta in the temperature case are linearly suppressed. 
In order to investigate the variation of the self-energy with the temperature near the 
critical temperature T c , we adopt the same method as in the case of finite density 
in the previous section. We only sketch the calculation. The critical temperature 
and the momentum dependence of the self-energy is obtained from the linearized 
equation, i.e. 



t 



k 2 \ 1 

V(k,p)f T (k 2 ), (28) 



V 



2T r 



2yk 2 



Analogous to the density case, we approximate the self-energy for temperatures close 
to T c by 

V(f) w S (1T) (T C ) f T {f) , (29) 
and obtain from the to next to leading order linearized Dyson- Schwinger equation 



S 



J (1T) 



d 3 k 



-k/T c 



X 



X 




(2tt) 3 (1 + e - fc / T 2 
d 3 k T r - T r e~ 2k ^ - 2ke- k ^ 



1/2 



k 3 (l + e- fc / T 2 



T^k^DZik 2 )^ 



(30) 



-1/2 



where := yk 2 . We find a Gaussian critical exponent for the chiral phase transition 
at finite temperature. 

5. Comparison of the phase transition at finite temperature and finite density 
All momentum integrals which occur in the previous sections, are ultra-violet finite, 
which is merely a consequence of asymptotic freedom and therefore independent of 
the details of the models under investigations. However, we have tacitly assumed 
that the integrals are also infra-red finite. In the case of finite density (section 3), 
this is of course true due to Pauli blocking of momentum states below the Fermi- 
momentum kf. In contrast, the case of finite temperature (section 4) needs further 
studies. 

Expanding the term in the second line of (|30[) with respect to small k, i.e. 



(31) 



T c - T c e~ 2k ^ - 2ke~ k ^ 1 2 

k 3( 1 + e -k/T c y ~ 12T 2 + °( k 



one observes that both integrals in ([30]) are infra-red finite, provided the product of 
vertex function and gluon propagator Ti(0, k 2 )D°^{k 2 ) does not diverge faster than 
1/k 2 . In the derivation of (0), we must exclude models, which incorporate quark 
confinement by the gluon propagator via infra-red slavery (see (|j)). For the gluon- 
propagator (|J), the infra-red divergence is screened by a tiny, but non- vanishing 
self-energy at temperatures somewhat below the critical one. This indicates that 
the restoration of the chiral symmetry is described by a different law than the simple 
one with critical exponent 1/2 in (p0[). Since the restoration of the chiral symmetry 
at finite density is completely insensitive to the low energy behavior of the interaction 
due to Pauli's principle, we expect a completely different behavior of the self-energy 
at finite density and at finite temperature, respectively, close to the critical point, 
if the model's interaction is infra-red sensitive. In contrast, we have shown that 
chiral transition at finite density is of the same nature as the transition at finite 
temperature for interactions which are "well behaved" in the infra-red region. 
Finally, we search for a relation between the critical Fermi-momentum and the crit- 
ical temperature. To this aim, we compare equations (|19|) and (|28|), which provide 
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the critical Fermi momentum and the critical temperature, respectively. Since in the 
latter case the tanh k/2T c rapidly approaches 1 for k 3> T c , and since the integral 
kernel in both equation strongly peaks at k « p, one finds that the function f{p 2 ) 
(103) coincides with /t(p^) ( P%D for p 2 S> T c 2 . On the other hand, both functions are 
normalized to 1 at p 2 = 0. We might therefore approximate 



(32) 



f(f) « /t(T) 

in the entire momentum range. Equating ( |i9|) and (|28| ) at p 2 = then yields 

roo roo { k \ 

J k dkkT^k^D^k 2 )/^ 2 ) = j o dk k tanh ( — J T x (0, k 2 )D™(k 2 )f(k 2 ). 

(33) 

A relation between k c and T c can be obtained, if one assumes that the function 
fcTi(0, k 2 )D a ^(k 2 )f(k 2 ) is monotonic decreasing with k. In this case, one finds that 



roo roo 

/ dkkT^k^DH^fXk 2 ) > / dkkT^k^D^k 2 )/^ 2 ) 

J kc J Icq 



where kn must be calculated from 



dk 



tanh 



2T r 



9{k - k ) 







with 8(k) the step-function. From this equation, we obtain ko 
fore the final result 

k r < 2 ln2T r w 1.386 T r 



(34) 



(35) 



2 In 2 T r and there- 



(36) 



Its instructive to confront this inequality with lattice data. Recent lattice simu- 
lations |T8| yield T c w 260 MeV. From (06|), we find A; r w 360 MeV. This value is 



in good agreement with the standard estimate pOf A; c w mjv/3, where is the 
nucleon mass. 

5. Conclusions 

We have studied a wide class of effective one-gluon exchange models. The basic 
assumption is a separation of the fermionic and gluonic energy scales justifying 
an instantaneous approximation. The quark interaction is constrained in order to 
reproduce the well-know behavior of the perturbative Green's functions at high 
momentum transfer. With some further assumption, which we believe to be justified 
for most of the one-gluon exchange models, we have been able to extract quite 
detailed relations for the critical behavior of various quantities, like the quark self 
energy S, at the chiral phase transition. We find that S as function of the Fermi 
momentum kf behave like 



S(p 2 = 0) oc 



\ 



k 2 

1 L 

kl 



(37) 
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where k c is the critical value of kf, and that the restoration of the chiral symmetry 
is insensitive to the infra-red behavior of the interaction due to Pauli blocking. In 
the case of finite temperature, the analogous behavior, i.e. 



is found, only if the interaction diverges weaker than 1/k 2 at small momentum 
transfer k. In the opposite case of an "infra-red sensitive" interaction, a completely 
different behavior of the self-energy near the critical density and the critical tem- 
perature, respectively, is expected. In the former case, we obtain Gaussian critical 
exponents supporting the recent ideas of Kocic and Kogut Q . Assuming further that 
the effective quark interaction is a monotonic decreasing function of the momentum 
transfer, we obtain the inequality 



which is in good agreement with recent QCD lattice simulations. 
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